SHARP L 1 STABILITY ESTIMATES FOR 
HYPERBOLIC CONSERVATION LAWS 



Paola Goatin 1 ' 2 and Philippe G. LeFloch 1 

ABSTRACT. We introduce a generalization of Liu- Yang's weighted norm to linear and to nonlinear hy- 
perbolic equations. Following an approach due to the second author for piecewise constant solutions, 
we establish sharp L 1 continuous dependence estimates for general solutions of bounded variation. Two 
different strategies are successfully investigated. On one hand, we justify passing to the limit in an L 1 
estimate valid for piecewise constant wave-front tracking approximations. On the other hand, we use the 
technique of generalized characteristics and, following closely an approach by Dafermos, we derive the 
sharp L 1 estimate directly from the equation. 



1. Introduction 

We are interested in the continuous dependence of entropy solutions to hyperbolic conservation 
laws 

d t u + d x f(u) = 0, u{x,t) £R, x eR, t > 0, (1.1) 

where the flux / : R — > R is a smooth and convex function. After works by Liu and Yang [22] and 
Dafermos [9], we aim at deriving sharp L 1 estimates of the form 



ii i 



(t) - u 1 {t)\\ w(t) + / M(T;u I ,u II )dT<\\u n (s)-u I (s)\\ wia) , 0< s <t, (1.2) 



I u 

for any two entropy solutions of bounded variation u 1 and u 11 of (1.1), where ||-|| TO (t) is a weighted 
L 1 -norm equivalent to the standard L 1 norm on the real line. In (1.2), the positive term M(t; u 1 \u TI ) 
is intended to provide a sharp bound on the strict decrease of the L 1 norm. The estimate with w = 1 
and M = is of course well-known. 

Recall that the fundamental issue of the uniqueness and continuous dependence for hyperbolic 
systems of conservation laws was initiated by Bressan and his collaborators (see [1, 2] and the 
references therein). A major contribution came from Liu and Yang [22, 23] who introduced a 
decreasing L 1 functional ensuring (1.2) for scalar conservation laws and systems of two equations. 
This research culminated in papers published simultaneously by Bressan, Liu and Yang [4], Hu and 
LeFloch [14], and Liu and Yang [24], which contain particularly simple proofs of the continuous 
dependence of entropy solutions for systems. 

In the present paper, we restrict attention to scalar conservation laws and, following the approach 
developed by the second author (Hu and LeFloch [14] and LeFloch [18]), we investigate the stability 
issue from the standpoint of Holmgren's and Haar's methods ([21] and the references therein). The 
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problem under consideration is (essentially) equivalent to showing the uniqueness and L 1 stability 
for the following hyperbolic equation with discontinuous coefficient: 

d t i/) + d x (ai>) =0, il>(x,t) eR, x e R, t > 0. (1.3) 

That is, for solutions with bounded variation we aim at deriving an estimate like 

\\m\\w(t)+ f M(t; a, ^)dr< W(s)\\ wW , 0<s<t. (1.4) 

J s 

For the application to (1.1) one should define a by 

/ / /In fiu 11 )- f{u r ) 

a = a(u\u H ) = — — 4 — — T — '-. (1.5) 
u 11 — u 1 

One may also consider the equation (1.3) for more general coefficients a. 

Recall that the existence and uniqueness of solutions to the Cauchy problem associated with (1.3) 
was established in LeFloch [16] in the class of bounded measures, under the assumption a x < E for 
some constant E. The latter holds when a is given by (1.5) (at least when u 1 and u 11 contain no 
rarefaction center on the line t = which holds "generically" ) . See also Crasta and LeFloch [5, 6] 
for further existence results. 

It must be observed that we restrict attention here to more regular solutions, having bounded 
total variation, as this is natural in view of the application to the conservation law (1.1). In this 
direction, recall that an L 1 stability result like (1.4) was established in [14] (see, therein, Section 5, 
and our Theorem 2.2 below) in the class of piecewise Lipschitz continuous solutions, with M = 
however. This uniqueness and stability result was achieved under the assumption that the coefficient 
a does not contain any rarefaction shock (see Section 2 below for the definition). In [14], the following 
essential observation was made: 

The linearized equation (1.3)-(1.5) based on two entropy 

solutions of (1.1) does not exhibit rarefaction shocks. ^ ^ 

(This is also true for systems of conservation laws, as far as solutions with small amplitude are 
concerned.) One of our aims here is to extend the L 1 stability result for (1.3) in [14] to arbitrary 
solutions of bounded variation. 

The present paper relies also heavily on the contribution by Liu and Yang [22] who, for ap- 
proximate solutions constructed by the Glimm scheme, discovered a weighted norm having a sharp 
decay of the form (1.2). Subsequently, the Liu- Yang's functional was extended by Dafermos ([9], 
Chapter 11) to arbitrary functions of bounded variation (BV) and, using the notion of generalized 
characteristics, Dafermos derived precisely an estimate of the form (1.2) valid for BV solutions. 

The aim of this paper is to provide a new derivation and some generalization of this L 1 functional. 
Toward the derivation of bounds like (1.2) or (1.4), we make the following preliminary observations: 

(1) The geometrical properties of the propagating discontinuities in a (Lax, fast or slow under- 
compressive, rarefaction shock, according to the terminology in [14]) play an essential role. 
It turns out that the (jump of the) weight w(x, t) should be assigned precisely on each under- 
compressive discontinuity. On the other hand, Lax discontinuities are very stable and do not 
require weight, while (in exact entropy soutions) rarefaction shocks do not arise, according 
to (1.6). 

(2) Certain (invariance) properties on the coefficient a are necessary to define the weight globally 
in space; see (2.9)-(2.10) in Section 2. 

(3) The weight however is far from being unique and we believe that this flexibility in choosing 
the weight may be helpful in certain applications. 
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The content of this paper is as follows. 

In Section 2, we consider piecewise constant solutions of (1.3) and introduce a class of weighted 
norms satisfying a sharp bound of the form (1.4). See Theorem 2.3 below. All undercompressive 
and Lax discontinuities contribute to the decrease of the L 1 norm. For the sake of comparison, we 
also consider the L 1 norm without weight; see Theorem 2.2. 

In Section 3, we point out that the setting of Section 2 covers the case of the conservation law 
(1.1). Passing to the limit in wave front tracking approximations, in Theorem 3.5 we arrive to 
the sharp bound (1.2) for general BV solutions. The proof is based on fine convergence properties 
established earlier by Bressan and LeFloch [3] and on a technique of stability of nonconservative 
products developed by DalMaso, LeFloch, and Murat [10] and LeFloch and Liu [20]. 

Next, in Sections 4 and 5 we return to the equation (1.3) studied in Section 2 but, now, we deal 
with general BV solutions. We follow closely ideas developed by Dafermos [8, 9] for solutions of 
(1.1), and extend them to the linear equation (1.3). Using generalized characteristics we establish 
first a maximum principle in Theorem 4.5. Finally, in Theorem 5.1 using the technique of generalized 
characteristics, we establish the sharp L 1 stability property (1.4) directly, for general BV solutions 
of (1.3). The result applies in particular to the conservation law (1.1) and allows us to recover (1.2). 

Throughout the paper, we always assume that all functions of bounded variation under consider- 
ation are normalized to be defined everywhere as right-continuous function. 

2. Decreasing Weighted Norms for Piecewise Constant Solutions 

Given a piecewise constant function a : R xR + — » R, let us consider the linear hyperbolic equation 

d t ^ + d x (aip) =0, iP(x,t)eR, (2.1) 

and restrict attention to piecewise constant solutions. By definition, the function a admits a set of 
jump points J (a), consisting of finitely many straightlines defined on open time intervals, together 
with a finite set of interaction points 1(a), consisting of the end points of the lines in J (a). The 
function a is constant in each connected component of the complement C(a) of 2(a) U J (a). At a 
point (x,t) € J (a) we denote by A a = X a (x,t) the speed of the discontinuity and a± = a±(x,t) = 
o(x±, t) the left- and right-hand traces. It is tacitly assumed that the discontinuity speeds A a remain 
uniformly bounded. Finally the function is normalized to be right-continuous. A similar notation is 
used for the function tp. 

The geometrical property of the coefficient a play a central role for the analysis of (2.1), so we 
recall the following terminology [14]: 

Definition 2.1. A point (x,t) G J(a) is called a Lax discontinuity iff 

a-(x,t) > X a (x,t) > a + (x,t), 
a slow undercompressive discontinuity iff 

\ a (x,t) < mm(a-(x,t),a + (x,t)), 
a fast undercompressive discontinuity iff 

\ a (x,t) > m&x.(a-(x,t), a + (x,t)) , 
and a rarefaction- shock discontinuity iff 

a-(x,t) < X a (x,t) < a + (x,t). 

For each t > 0, we denote by C(a),S(a),J-(a), and 7Z(a) the set of points (x, t) G J(a) correspond- 
ing to Lax, slow undercompressive, fast undercompressive, and rarefaction-shock discontinuities, 
respectively. 
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Theorem 2.2. Consider a piecewise constant speed a = a(x,t). Let ip be any piecewise constant 
solution of (2.1). Then we have for all < s < t 



\m)hHR)+ f E 2(a-(x,T)-\ a (x,T))\i>-(x,T)\dT 
= U(s)\\ LHR) + f Yl 2(X a (x,T)-a.(x,T))\^(x,r)\dT. 



(x,T)e1Z(a) 



In (2.2), the left-hand traces are chosen for definiteness only. Indeed it will be noticed in the 
proof below that for all (x, r) <G C(a) U lZ(a) 

(A Q (x,r) - a_(x,r)) \iP-(x,r)\ = -(X a (x,r) - a+(x,r)) \^+(x,r)\ 

Observe that the Lax discontinuities contribute to the decrease of the L l norm, while the rarefaction- 
shocks increase it. On the other hand, the undercompressive discontinuities don't modify the L 1 
norm. When a contains no rarefaction shocks (this is the case when (2.1) is a linearized equation 
derived from entropy solutions of a conservation law, as discovered in Hu and LeFloch [14]), Theorem 
5.1 yields 

\m)hHR) < \Ms)\\ LHR) , 0<s<t, (2.3) 

where we neglected the favorable contribution of the Lax discontinuities appearing in the left-hand 
side of (2.2). In particular, (2.3) implies that the Cauchy problem for (2.1) admits a unique solution 
(in the class of piecewise constant functions at this stage), provided a has no rarefaction-shock 
discontinuities. 

On the other hand, it is clear that the sign of the function ip is important for the sake of deriving 
the L 1 stability of the solutions tp of (2.1). For instance, if -0 has a constant sign for all (x,t), then 
(2.3) holds as an equality 

\m)\\LHR) = W(8)\\ L HR), 0<s<t, 

which implies that the Cauchy problem for (2.1) admits at most one solution ip of a given sign. 

Proof. Denote by P(E) the projection of a subset E of the (x,t)-plane on the t-axis. By definition, 
any piecewise Lipschitz continuous solution ip is also Lipschitz continuous in time with values in 
^(R). So, it is enough to derive (2.2) for all t g E := P(X(a) Ul(ip)). The latter is just a finite 
set. The following is valid in each open interval / such that I D E = $. 

We denote by Xj(t) for t £ I and j = 1, ■ ■ • ,m the discontinuity lines where the function ip(.,t) 
changes sign, with the convention that 

ij)(x,t) > for x € [xj(t),x j+ i(t)]. (2.4) 



Set ipf(t) = ip±(xj(t),t), Xj(t) = X a (xj(t),t), etc. Then by using that ip solves (2.1) we find (for all 
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t in the interval /) 
It 



x, t) \ dx 



R 



= ^E(" 1 ) J / rl>(x,t)dx 

= E(- 1 ) j ( r +l(t) dMx,t)dx + X j+1 (t)^r +1 ( t ) - X^^it) ) 

j=l \ J Xj{t) J 

m / » x (t) ^ 

= E(- 1 ) j / -d x {a{x,t)i>{x,t))dx + A i+ i(t)V7 +1 (t) - \j(t)tf(t) 

j=i yJxjit) j 

m 



The Rankine-Hugoniot relation associated with (2.1) reads 

(a+(t) - Xj(t))^(t) = (aj(t) - X,(t))^-(t), (2.5) 

therefore by (2.4) 

7 n m 

- J^(x,t)\dx = 2j2M4(t)-^(t))\^f(t)\. (2.6) 

Consider each point Xj(t) successively. If Xj(t) is a Lax discontinuity, then a~(t) > Xj(t) > a+(t) 
and both coefficients ±(a • (t) — Aj(t)) are negative. If Xj(t) is a rarefaction-shock discontinuity, then 
a~(t) < Xj(t) < a + (t) and the coefficients ±(a^(t) — Xj(t)) are positive. These two cases lead us 
to the two sums in (2.2). Indeed one just needs to observe the following: if (x,t) correspond to a 
Lax or rarefaction-shock discontinuity of the speed a, but ip does not change sign at (x, r) (so it is 
not counted in (2.6)), then actually by the Rankine-Hugoniot relation (see (2.5)) we conclude easily 
that 

ip-{x,T) = ip + (x,r)) = 0, 

and so it does not matter to include the point (x,r) in the sums (2.2). 

Suppose next that Xj(t) is an undercompressive discontinuity. Then the two sides of (2.5) have 
different sign, therefore 

(at(t) - XjiW+it) = (aj(t) - \j(t))il>j(t) = 0, 
and the corresponding term in (2.6) vanishes. □ 

Our objective now is to derive an improved version of Theorem 2.2, based on a weighted L 1 norm 
adapted to the equation (2.1). For piecewise constant functions, we set 

IIV'WIU) : = / \4>(x,t)\w(x,t)dx, (2.7) 
Jr 

where w = w(x,t) > is a piecewise constant and uniformly bounded function. We determine this 
function based on the following constrain on its jumps, at each discontinuity of the speed a, 

' < if (x,t) G S(a), 

(2.8) 



w+(x, t) — W- (x, t) 



> if (x,t) € T{a). 
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The weight is chosen so that the left-hand trace of a slow undercompressive discontinuity and the 
right-hand trace of a fast one are weighted more. This is consistent with the immediate observation 
that the terms (\j(t) — aj (t)) \tpj(t)\ and (a^(t) — Xj(t)) \ip~j~(t)\ have a favorable (negative) sign 
for slow and fast undercompressive discontinuities, respectively. On the other hand, the jumps of w 
at Lax or rarefaction-shock discontinuities will remain unconstrained. This choice is motivated by 
the two observations: (i) Lax shocks already provide us with a good contribution in (2.2), and (ii) 
rarefaction shocks are the source of instability and non-uniqueness and cannot be "fixed up" . 

The constrain in (2.8) is different for slow and for fast undercompressive discontinuities. To 
actually exhibit a (uniformly bounded) weight satisfying (2.8), we put a restriction on how the 
nature of the discontinuities change in time as wave interactions take place. (An incoming wave 
may be a slow undercompressive one and become a fast one after the interaction, etc. A different 
constrain is placed before and after the interaction.) 

Precisely, we suppose that, to the speed a = a(x,t), we can associate on one hand a function 
k : R x R + — > R having bounded total variation and such that J(n) C J {a) and T(n) C 2(a), and 
on the other hand a partition of the discontinuities 



so that, for each (x,t) G J (a), the limits n± = K±(x,t) determine if the wave is slow or fast on its 
left or right side, as follows: 



Here we use sgn (y) = —1,0,1 iff y < 0,y = 0,y > 0, respectively. Therefore a discontinuity 
(x,t) 6 J 1 {a) (for instance) is 



Furthermore, to measure the strength of the jumps, we introduce a piecewise constant function, 
b = b(x, t), having the same jump points as the function a. For instance, we could assume that there 
exist constants C\ , C2 > such that at each discontinuity of a 



J (a) = J 1 (a) U J 11 (a), 



(2.9) 




(2.10) 



a Lax one iff k_ < and k + > 0, 
a slow undercompressive one iff K- > and n + > 0, 
a fast undercompressive one iff k_ < and k + < 0, 
a rarefaction-shock iff k_ > and n + < 0. 



C 1 \a + (y,t)-a.(y,t)\ < \b+(y,t) - b_(y,t)\ < C 2 \a + (y,t) - a.(y,t)\. 



(2.11) 



However, strictly speaking, this condition will not be used, in the present section at least. 
Based on the functions k and b and for t except wave interaction times, we can set 



V I (x,t)= Yl \b+(v,t)-b-(y,t)\, 



(y,t)e^ 7 (a), 



(2.12) 



V n (x,t)= £ \b+(y,t)-b-(y,t)\, 



(s,()eJ"(»), 

y<x 



so that the total variation of b(t) on the interval (— 00, x) decomposes into 



TV* 00 (b(t)) = V I (x,t) + V II (x,t). 



(2.13) 
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Fix some parameter m > 0. Consider now the weight-function defined for each (x, t) G C(a) by 

{m + y 7 (oo, i) - F 7 (x, t) + V u (x, t) if k(x, t) > 0, 
(2.14) 
m + V 7 (x, t) + V n (oo, t) - V n (x, t) if k(x, t) < 0. 

It is immediate to see that indeed (2.8) holds and that with (2.11) 

m < w(x,t) < m + TV{b{t)) < m + C 2 TV(a(t)), x G R. (2.15) 

Note also that the weight depends on b and a, but not on the solution. 

Theorem 2.3. Consider a piecewise constant speed a = a(x,t) admitting a decomposition (2.9)- 
(2.10) and satisfying the total variation estimate (2.15). Consider the weight function w = w(x,t) 
defined by (2.13). Let tp be any piecewise constant solution of the linear hyperbolic equation (2.1). 
Then the weighted norm (2.7) satisfies for all < s <t 

ll^)IU(t) 

+ /E (2m + TV(6)-|6 + (x,r)-6_(x,r)|)|a_(ar,r)-A(x,r)|^_(x,r)|dr 
+ 1 \ b +(x,r) -6_(x,r)| |a_(x,r) - X(x,t)\ \ip-(x,r)\dT 

Js (x,r)65(a)Uf(«) (2-16) 



Ms)\Us)+ [ E (2m + TV(b))\a.(x,r) - X(x,T)\\^(x,r)\dT 
+ E \ b +( x , T ) -b-(x,T)\\a-(x,T) - \(x,T)\\ip-(x,T)\dT. 

Js / \^T-,r-.\ 



(x,r)en(a) 

The statement (2.16) is sharper than (2.2), as all discontinuities contribute now to the decrease 
of the weighted L 1 norm. Note that asm^oo, we recover exactly (2.2) from (2.16). 

Proof. We proceed similarly as in the proof of Theorem 2.2. However, xj(t) for t E I (some open 
interval avoiding the interaction points in a or tp) denote now all the jump points in either a or tp. 
We obtain as before the identity 

d_ 

di 



ip(x, t) | w(x, t) dx 

= £ ((\ 3 (t) - a-(t)) \^(t)\w-(t) + (a+(t) - A,-(t)) \^+(t)\w+(t)) ^ 

.7 = 1 
777 

= £ (sgn {Xj(t) - aj(t)) w~(t) + sgn (a+(t) - \j(t)) w+(t)) \Xj(t) - aj(t)\ ^ (t)\ , 

j=i 

where we used the Rankine-Hugoniot relation (2.5). 

If xj(t) is a Lax discontinuity in J 1 (a), then by (2.11) we have k_ < and k + > 0. So by (2.14) 
we find 

vj- =m + V'ixjit)-) + V n {oo) - V H ( Xj (t)-), 
w + = m + V T (oo) - V I (x j (t)+) + V II (x j (t)+), 
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and so 

sgn (Xj(t) - aj(t)) w~(t) + sgn (a+(t) - \j(t)) w+(t) 

= -wj{t) -w+(t) (2.18a) 
= -2m-TV(b) + \b+(t)-bj(t)\. 

If Xj(t) is a rarefaction-shock discontinuity in J 1 (a), then by (2.11) we have K- > and k + < 0. 
By (2.13) we find 

wj =m + V\oo) - V' (£,(*)-) + V n ( Xj (t)-), 
w + = m + V I (x j (t)+) + V H (oo) - V II (x j (t)+), 

and so 

sgn (Xjit) - aj(t)) wj(t) + sgn (a+(t) - Xj(t)) w+(t) 

= w -(t)+w+(t) (2.18b) 
= 2m + TV(b) + \bf(t)-bJ(t)\. 

If Xj(t) is a fast undercompressive discontinuity in J 1 (a), then by (2.11) we have k_ < and 
k + < 0. By (2.13) we find 

sgn (\j{t) - aj(t)) wj(t) + sgn (a+(t) - Xj(t)) w+(t) 
= w-(t)-w+(t) 

= m + V^Xjit)-) + V H {oo) - V H (xj(t)-) (2.18c) 

- m - V I (x j (t)+) - V H (oo) + V II {x j {t)+) 
= -\b+(t)-b-(t)\. 

Similarly for slow undercompressive discontinuities in J 1 (a) we obtain 

sgn (Aj (t) — aj (t)) wj (t) + sgn ( a + (t) - X,(t)) w+ (t) = -\b+ (t) - bj (t)\. (2.18d) 

Using (2.18) in (2.17) we conclude that 



\\m\Ut) 




which is equivalent to (2.16). □ 

Using that IZ(a) is included in the set of points where ip changes sign, it is easy to deduce from 
(2.16) that: 
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Corollary 2.4. Under the assumptions and notations in Theorem 2.3, we have for all < s < t 

\\m\Ut) 

<\M8)\\ wW +(2m + TV(bj) sup \b + (x,r)-b-(x,r)\ f TV(^(r))dr 

S<T<t 

and, in particular, letting m — > oo 

\Mt)\\LHR)<\Ms)\\ LHR) +2 sup \b + (x,T)-b.(x,r)\ [ TV^(r))dr. (2.19) 

S<T<t 

Finally, in view of Corollary 2.4, in case the function a contains no rarefaction shocks, we deduce 
that 

\m)\Ut) < u(s)\\ w{s) , o< s <t. 

Observe that this result is achieved, based on a weight that depends on an arbitrary function, 
b, and on the sole assumption that a decomposition (2.9)-(2.10) of the jumps of a is available. 
However, our result in this section covers only piecewise constant solutions. We will see in Section 
5 that a stronger structure assumption on the coefficients a is necessary to handle general solutions 
of bounded variation. 

3. Sharp L 1 Estimate for Hyperbolic Conservation Laws 

In this section, we apply Theorem 2.3 to the case that a is the averaging coefficient (1.5) based 
on two entropy solutions of (1.1). First, we check that the assumptions required in Section 2 on 
the coefficient a do hold in this situation. Therefore Theorem 2.3 applies to the piecewise constant 
solutions defined by the wave-front traking (also called polygonal approximation) algorithm proposed 
by Dafermos in [7]. Next, we observe that, with a suitable choice of the definition of the wave 
strengths, the weighted norm in Section 2 reduces to Liu- Yang's functional. Finally we rigorously 
justify the passage to the limit in the estimate of Theorem 2.3 when the number of wave fronts tends 
to infinity and exact entropy solutions of (1.1) are recovered. 

Consider the nonlinear scalar conservation law: 

8 t u + d x f(u) = 0, u(x,t)eR, (3.1) 

where the flux / : R — > R is a smooth function. Let u 1 and u 11 be two bounded entropy solutions of 
(3.1) having bounded total variation. Given h > let us approximate the data u 1 (0) and u 11 (0) by 
piecewise constant functions u I,h (0), u II,h (0), having finitely many jumps and such that as h — > 

^(0)^^(0), u n ' h (0) -» u u (0) in the L 1 norm, (3.2) 

TVV^O)) -f TViu 1 '(0)), TV(u n ' h (0)) -» TV{u H (0)). (3.3) 

Applying Dafermos' scheme [7], we can construct corresponding, piecewise constant, approximate 
solutions u 1 ' 11 and u II,h having finitely many jump lines and for t > s > and p£ [1, oo] 

ll^ h (t)|| LP(jR ) < ||^(s)|k,0R), \\u n > h (t)\\ LP(R) < \\u n ' h (s)\\ LPm , (3.4) 
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and for all -oo < A + M (t - s) < B - M (t - s) 

TV^^{u^{t)) < TV ^ (u I,h (s)), 
TVEZg-${u"*(t)) < TV*(u n ' h (s)). 

More precisely, the functions u I,h and u II,h are exact solutions of (3.1) satisfying therefore the 
Rankine-Hugoniot relation at every jump. They contain two kinds of jump discontinuities: Lax 
shocks satisfy the so-called Oleinik entropy inequalities, while rarefaction jumps do not, but have 
small strength, that is 

\u I ' h (x+,t)-u I ' h (x-,t)\<h, \u H ' h (x+,t)-u n ' h (x-,t)\<h. (3.6) 

Furthermore, for a subsequence h — > at least, we have for each time t > 

t/' h (t) -f u J (i), u n ' h {t) -» u n {t) in the L 1 norm. 

To study the L 1 distance between these approximate solutions, we set 

^ : = u n > h - u 1 ^, 

which is one solution of the linear hyperbolic equation 

First of all, based on Theorem 2.2 and (3.5)-(3.6), we obtain immediately: 

Theorem 3.1. The approximate solutions u I)h and u H ' h satisfy the following L 1 stability estimate 
for all0<s<t 

\\u n ' h (t)-u^ h (t)\\ L , {R) 



+ Yl 2(a h {x-,T) - X ah {x,T))\u n ' h (x-,T) -^^{x-^^dr 

Js {x,r)€C{a) 

(s) - ^ h ( S )|| L1(jR) +2h(t- s) II/" |U (TV(u^ h (0)) + TV(u n > h (0))) . 



< \\u H > h 

From the functions u 1 and u 11 we define the function a as in (3.7). Recall that the wave front 
tracking scheme converge locally uniformly (see the proof of Theorem 3.5 below for a the definition), 
so that the BV solutions u 1 and u 11 are endowed with additional regularity properties. Consider for 
instance the function u 1 . In particular, for all but countably many times t and for each x, either x 
is a point of continuity of u 1 in the classical sense (say (x,t) G C(u r )) or else it is a point of jump 
in the classical sense (say (x,t) G J(u r )) and, to the discontinuity, one can also associate a shock 
speed, denoted by A 7 (x,t). 

From the properties shared by u 1 and u 11 , one deduces immediately a similar property for the 
coefficient a. Excluding countably many times at most, at each point of jump of a we can define the 
propagation speed X a (x,t) of the discontinuity located at the point (x,t). Namely, we have 



X a (x,t) 



A 7 (x,t) (x,t) £ Jiu 1 ), 
X n (x,t) if (x,t) G J(u H ). 

In the limit ft^Owe deduce from (3.8) that: 
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Corollary 3.2. For all < s < t we have 
||u"(t)-</(t)||z,H*) 

+ [ ^(a(x-,T)-X a (x,T))\u II (x-,T)-u I (x-,T)\dr (3.9) 

Js (x,r)eC(a) 
< \\u II (s)-u I (s)\\ Llm . 

We omit the proof of Corollary 3.2 as (3.9) is a consequence of a stronger estimate proven in 
Theorem 3.5 below (by taking m — ► 00 in (3.15)). Note that (3.9) is a stronger statement than the 
standard L 1 contraction estimate 

\\u H (t) - </(i)|| L l (jR) < \\u H ( S ) - u J ( S )|| L1(J?) . 

Proof. We apply the estimate (2.2) with tp replaced with u 1I,h — u I,h . We just need to observe (see 
[14]) that all the rarefaction-shock discontinuities in a h are due to rarefaction fronts in u I,h or in 
u II,h , which have small strength according to (3.6). In other words we have 



f E 2(A a (x,r)-a_(x,r))|V-(x,r)|dr 

< sup 2 \a + (x,r) - a_(x,r)| f TV(iP(t))oIt 

S<T<t 

<2\\f"\\ 00 h J*TV(^(r))dT 

< 2 ||nU h(t-s) (TViu^iO)) + TV(u n ' h (0))) . 

This establishes (3.8). □ 

We now want to apply Theorem 2.3 and control a weighted norm of u 1I,h — u I,h . In this direction, 
our main observation is: 

Lemma 3.3. When the function f is strictly convex, the coefficient a h satisfies all of the assump- 
tions (2.9)-(2.10). 

Proof. The function a h is piecewise constant, and we can associate to this function an obvious 
decomposition of the form (2.9). To establish (2.10), consider for instance a jump point (x,t) £ 
J{u I,h ) n C(u II,h ), together with its left- and right-hand traces uL and u+. Since u I,h is a solution 
of (3.1), the corresponding speed A = X(x,t) satisfies the Rankine-Hugoniot relation: 

-\{u I + -u I _)+f(u I + )-f(u I _) = 0. 

Thus the term in the left-hand side of (2.10) takes the form 

/(«") - /(i4) /«) - f{ul) 



a±(x, t) — X(x, t) 



i 11 — u ± u + ~ u - 



£ (f'{0u n + (1-9) 4) - f'{0v^ + (1-0) 4)) d6. 
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Thus we obtain 

a±(x,t) — X(x,t) = n (u 11 — u^, 

f 1 f 1 / \ (3.10) 

M: =y J f"(p(0u n + (1-0)4) + (1- p)(eu I T + (l-6)u I ± ))ed0dp. 

Since / is strictly convex, the coefficient is bounded away from zero. In view of (3.10), if we now 
choose K,(x,t) := u II,h — u T ' h , the desired property (2.10) holds true. □ 

Next, we define the weight w h = w h (x,t) associated with the function a , by the formula (2.14) 
in which we specify 

K h (x,t):=u II ' h -u I ' h . (3.11) 
It follows immediately from Theorem 2.3 that: 

Theorem 3.4. Suppose that the function f is strictly convex. The approximate solutions constructed 
by Dafermos scheme satisfy the L 1 stability estimate for all < s <t 

h n,h {t) _ u l,h mwHt) 

+ 1 £ (2m + TV(b h )-\b h (x+,T)-b h (x-,T)\) 

S (x,T)eC(a h ) 

\a h (x-,T) - X h (x,T)\ \ U II ' h ( y X-,T)-U I ' h (x-,T)\dT 

+ f E \b h (x+,r) -b h (x-,r)\\a h (x-,r) - X h (x,r)\ 

S (x,T)eS{a h )\jF(a h ) (3-12) 
\u"' h {x-,T)-U^ h {x-,T)\dT 

= \\u II > h (s)-u I > h (s)\\ wHs) 

+ 1 hm + TV(b h ) + \b h (x+, T ) -b h (x-,T)\) 

\a h (x-,T) - X h (x,T)\ \ U II ' h (x-,T)-U I ' h (x-,T)\dT, 



(x,r)GTZ(a h ) 



where a h is the averaging coefficient defined in (3.7) and X h (x,T) represents the speed of the discon- 
tinuity located at (x,r) G J(a h ). 

We emphasize that (3.12) is an equality in which the contribution to the L 1 norm of each type of 
wave appears clearly. The coefficient a h exhibits three types of waves: the Lax and undercompressive 
discontinuities in a h contribute to the decay of the L 1 weighted distance. The statement (3.12) 
quantifies sharply this effect. On the other hand, the rarefaction-shocks appearing in the right-hand 
side of (3.12) increase the L 1 norm. 

In the rest of this section, we assume that the function b = b h is chosen to be specifically 

( u^ h (x+A) -u^ h (x-A) if (x,t) G J(u T ' h ), 
b h (x+,t)-b h (x-,t) = I rr ,\ ' \ TT \, ' ' , ttI\ (3-13) 

V ,J K ,J \ u H ' h (x+,t)-u H ' h (x-,t) if (x, t) G J(u ' ), K ' 

but a more general definition is possible. 

Our next purpose is to pass to the limit {h — > 0) in the statement established in Theorem 3.4 for 
piecewise constant approximate solutions. We recover here a result derived by Dafermos [9] via a 
different approach. Recall the notation C(u r ), S^u 1 ), etc introduced earlier. Denote by Tiu 1 ) the 
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countable set of interactions times. Let V 1 ^) be the total variation function associated with u T (t). 
Based on the functions V T (t) and V 11 ^), we then define the weight w as in (2.14) but with (2.12) 
replaced by the total variation functions of u 1 (t) and u 11 (t), with k := u 11 — u 1 and 



b(x + ,t)-b(x-,t) = { u n (x+it) _ u i I{x _ tt) i£{Xtt)e j iu ii ) . ^ 



Furthermore, to any functions of bounded variation u, v, w in the space variable x (the time 
variable being fixed) we associate the measure on R 

H = (a(u, v) — f'(u)) (v — u) dw 

understood as the nonconservative product in the sense of Dal Maso, LeFloch and Murat [10] and 
characterized by the following two conditions: 

(1) If B is a Borel set included in the set of continuity points of w 

fi(B)= I (a(u,v) - f'(u)) (v -u)dw, (3.15a) 

J B 

where the integral is defined in a classical sense; 

(2) If x is a point of jump of w, then 



MM) = j (i a ( u +' v +) -a(u-,u+)) (v + -u + ) 

+ (a(u_, v_) — a(u_, -u + )) (u_ — u_)^ \w + — 



(3.15b) 



with u± = u(x±), etc. 

Note that, if u = u 1 and v = u 11 , the two terms (a(u±,v±) — a(u-,u + )) (v± — u±) in fact coincide. 

Theorem 3.5. Let the function f be strictly convex and let u 1 and u 11 be two entropy solutions of 
bounded variation of the conservation law (1.1). For all < s < t we have 

\\u H (t) -u'MWvw 

+ ^2 q\a{u I (x-),u II (x-)) -a{u I (x+),u I (x-))\\u II (x)-u I (x)\dT 

+ f q\a(u I (x-),u II (x-)) -a(u II (x+),u II (x-))\\u II (x)-u I (x)\dT 

Js (x,r)€C(a)nJ{u") (3.16) 

+ / / (aiu^u 11 )- f'iu 1 ))^ 11 -u^dV 1 dr 

Js JR 

+ J s J {a(u I ,u II )-f(u II ))(u I -u II )dV II dr 
< \\u H (s) - n J (s) |U( S ). 

where q = q{r) = 2 m + TV {u 1 {t)) + TV {u 1 1 (t)) . 



14 



P. GOATIN AND P.G. LEFLOCH 



Observe that the terms in integrals in (3.16) globally contribute to the decrease of weighted norm, 
as is better seen rewriting the formula as follows (Vj and Vj 1 being the continuous parts of the 
measures V 1 and V 11 ): 



w{t) 



-I ( q ~ ' n + ~ u -\) \ a { u -i u +) -a(«+,«_)| lu 11 -u T \dr 

(a:,T)e£(a)n:7(W) 

- / £ (? - W+ - ui'l) \a(u z _, ul 1 ) - a^yj) | \u H -u z \dT 

(i,r)e£(a)nj(u") 

- I ^ |a(t/,?/ 7 ) -a(u+,u')| \u n -u r \ \u z + -u r _\dT 

(x,T)e(s(a)UT(a))nJ(u I ) 

+ /'* Y, ) " ««,«-) | lu 11 - u r \ lui 1 - uL'ldr 

(^T)e(«S(a)UF(a))n i 7(u- rI ) 

+ / / ,u n ) - f'iu 1 )^ 11 - u^dV* dr 

Js Jr 

+ [ [ laiu^u 11 )- fiu 11 )]^ 1 -u n \dV c n dr 

Js JR 
< \\U II (S)-U I (8) 



i 
s: 

i; 



\w(s)- 



(3.16') 



The following estimate is a direct consequence of the definition (3.15): 

Lemma 3.6. There exists a constant C > such that for all functions of bounded variation 
u,u,v,v,w defined on some interval [a, (5\ 



rP r/3 
/ (a(u, v) — f'{uj) (v — u) dw — I (a{u,v) — /'(it)) (v — u) dw 

J a J a 

<C (\\u- u\\ Lx{oiy/3) + \\v- v\\ Lao{ayf3) ) TV [a ^](w). 



(3.17) 



Proof of Theorem 3.5. 

Step 1 : Preliminaries. 

For each t > 0, the functions V I,h (t) and V n ' h {t) associated with the wave front tracking ap- 
proximations u 1 ' (t) and u II,h (t) are of uniformly bounded variation as h — > 0. The measures dV I,h 
and dV II,h are also Lipschitz continuous in time (with constant independent of h) for the weak 
convergence, except at interaction points. On the other hand, interaction times in the limiting so- 
lutions are at most countable. Therefore, extracting subsequences if necessary, the measures dV I,h 
and dV II,h converge to some limiting (non-negative) measures, say: 



dV^ h (t) -» dV'it) dV u > h (t) -» dV u {t). 
By lower semi-continuity, we have at each time t 

dV 1 {t) < dV 1 (t) , dV n {t) < dV n it) , 



(3.18) 



(3.19) 
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and, in particular, at each (x, t) 

V z {x, t) < V J {x, t), V H {x, t) < V H {x, t). (3.20a) 

F 7 (+oo,i) - V z (x,t) < F'(+oo,i) - V z (x,t), 

~~ _ _ (3 20b) 

V JI (+oc, t) - V J/ (x, t) < V 11 (+oo, i) - F /J (x, t). 

Based on the functions V 1 (t) and V 7/ (i), on the coefficient k := u 11 — u 1 and on the function in 
(3.14), we can define a weight denoted by w, along the same lines as in (2.14). We will show that 
the left-hand side of (3.16) is bounded above by 

Wu 11 (t) - u 1 (t)\\ w{t) 

+ f q\a{u I (x-),u II (x-)) -a(u I (x+),u I (x-))\\u II (x)-u I (x)\dT 

S (a;,T)e£(a)n:7(u') 

+ f q\a(u I (x-),u II (x-)) - a(u n (x+),u n (x-))\ \u n (x) -u J (x)|dr , 3 ^ 

S (x,T)€C(a)nJ(u Ir ) 

+ f I (a(uV)-/V)) (u 11 -u^dV^y^dr 

J s JR 

+ f f (a(uV) - f'iu 11 ))^ 1 -u n )dV n ( yi T)dT 

Js JR 

where q := 2m + V 1 (+oo) + V II (+oo), and that (3.21) coincides with the desired upper bound 
\\u n (s) — n / (s)|| ?i ,( s ). The former statement is postponed to Step 5 below and we focus now on the 
latter. 

Fix some t > s > and rewrite (3.12) in the equivalent form 

»«"•*(*) -« J ' fc (t)IU* ( t) 

+ f Yl (2m + TV(b h )^\a h (x-,T) - A h (x,r)| \u n ' h (x-,r) -u J ' h (x- r)| dr 

+ f Y, \b h (x+,T)-b h (x-,T)\(a h (x-,T)-\ h (x,T))(u II > h (x-,T)-u I ' h (x-,T))dr 

S {x,T)eJ'{a h ) 

+ f Y \b h (x+,r) -b h (x-,r)\ {a h (x-,r) -X h (x,r)) (t/' h (x-,r) - u h ^(x-,t)) dr 
= \\u II > h (s)-u I > h (s)\\ wHs) 

+ ! (2m + TV(b h )^j \a h {x-,r) - X h (x,r)\ \u n ' h (x-,T) - ^(x- r)| dr, 

S (x,T)&n(a h ) 

(3.22) 

or, with obvious notations, 

\\ u ii,h {t) _ u i^ m +Q h + Q h = \\ u ii,h( 8 ) -u^WWn Hs) +nl (3.23) 



As the maximum strength of rarefaction fronts in u 1 ' 11 and u H ' h vanishes with h (see (3.6)) and 
rarefaction shocks in a h arise only from these rarefaction fronts (see (1.6)), we have 



ft§ -» as h -> 0. (3.24) 
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On the other hand, we can always choose the (initial) approximations at time s in such a way that 

w(s) = w(s) (3.25) 

and 

lim \\u n ' h (s) - ^(s)|U (s) = ||n /7 ( S ) - ^( S )|U (S) . (3.26) 

It remains to prove that the limit of the left-hand side of (3.22) is exactly (3.21). This will be 
established in the following three steps. 

Step 2 : We will rely on the local uniform convergence of the front tracking approximations 
(see Bressan and LeFloch [3]). For all but countably many times r we have the following properties 
for u 1 (as well as for u 11 ): 

(1) For each point of jump z of u 1 there exists a sequence z h — > z such that for each e > there 
exists 5 > such that 

lu^fx) -u I (z+)\ + lu 1 (x) - u 1 (z+)\ < e for all x - z h G (0,5), 

T h T T T H ( 3 - 27a ) 

\u I ' h {x)-u I {z-)\ + \u I {x)-u I {z-)\<e for all x - z h G {-5,0) 
and (clearly) 

lu^x) - u I (z+)\ < e for all x - z G (0,5), 

T T ( 3 - 27b ) 

\u (x) - u (z-)\ < e for all x - z G (-5,0). 

(2) For each point of continuity z of u 1 and for each e > 0, there exists 5 > such that 

\u I ^ h (x)-u I {z)\ + \u I {x)-u I {z)\<e for all x-z G (-6,6). (3.28) 

We also recall from [3] that, for all but countably many times t, the atomic parts of the measures 
V 1 and V 11 coincide with the one of V 1 and V 11 , that is for each y G R 

V I (y+,t) - V\y-,t) = V I (y+,t) - V^y-^), 
V II (y+,t)-V II (y-,t) = V II (y+,t)-V II (y-,t). 

Following LeFloch and Liu [20] who established the weak stability of nonconservative products 
under local uniform convergence, we want to show that 



n h 2 (r) := / (a(u I '\y,T),u II ' h (y,r)) - /V' fc (i/,r))) (u n ' h (y, r) - u^ h {y, r)) dV T ' h (y) 

JR 

+ [ (a(u I ' h {y,T),u II > h {y,T)) - /V' h (y,r))) {u z > h (y, r) - u H ' h (y, r)) dV n > h (y) 

JR 

^ (a(n J (y, r),u n {y, r)) - f{u\y, r))) (u H (y, r) - u\y, r)) dV 1 \y) 
[ {a{n I (y,T),u II ( y ,T))-f\n II (y,T)))(n I (y,T)-u II (y,r))dV n (y). 

JR 



(3.30) 
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By Lebesgue dominated convergence theorem and since a uniform bound in r and h is available, it 
will follow from (3.29) that 

n h 2 = fn h 2 {T)dr-^ ft (a(«V) - /V))(«" -u I )dV I (y,T)dr 

J s Js Jr. (3 30 ,) 

+ / / {a{u\u II )-f'(u II )){u I -u II )dV II {y,T)dT. 



R 



Given e > 0, select finitely many (large) jumps in u 1 or u 11 , located at yi,y 2 , . . . y n , so that 

\u I (x+)-u I (x-)\ + \u II (x+)-u II (x-)\<e. (3.31) 



3=1,2,... ,n 



To each yj we associate the corresponding discontinuity point y^ in u I,h or u 11 ' . To simplify the 
presentation we will focus on the case where yj < yj < y 3 +\ < y^ +1 for all j. The other cases can 
be treated similarly. In view of the local convergence property (3.27)-(3.28) and by extracting a 
covering of the interval [yo,yn], we have also 



I ' h (x)-u I (x)\ + \u H > h (x) -u n {x)\ < 26, x £ (l/J )2 / i+ i) C (y^+i). (3.32) 



In view of (3.30) we can construct functions u\ and -u 7/ that are continuous everywhere except 
possibly at the points yj and such that the following conditions hold with u replaced by either u 1 



or u 11 : 



TV(u e ;R\{ yi ,... ,y n }) < C TV (u;R\ {y u . . . ,y n }) , 
<Ce, TV(u-u e ;R\{ yi ,... ,y n }) < C e, 



\u — u f 



(3.33) 



where C is independent of e. 
Consider the decompositions 



and 



/l t~b /i Til 

/ (a(tt J '\ u 7 ^) - /' (u 1 ' 11 )) (u n ' h - u T ' h ) dV z > h = V / ■+^2 
f (a(u' s u") - / V)) («" -u')dV' = Y,( ■ ■ ■ + E / 



j =i ) ■ / i •' i y ; } 

Here yft = yo = — oo and = y„+i = +oo. Thus in (3.30) we have to estimate 



n£(r) = / (a^V^) - /V ,h )) (u"' h - u^) dV J ' h 

- / {a^M 1 ) ~ fin 1 )) (u 11 - u 1 ) dV 1 ( 3 - 34 ) 
JR 

ih | rph 



T n + T ft 



with 



n ,. 

T* : = V / (a(u*'\ u"' fc ) - f{u^ h )) (u H ' h - u z > h ) dV z > h 
j= i 

n 

-E / (a(u J ,«")-/V))(«"-« 7 )^ 
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and 

n „ 

T 2 h := £ / (a(u J , - f'(u z ' h )) (u n ' h - u T ' h ) dV r ' h 

-Y, / (a(u^u")-/V))(«"-«')<*P'- 
First, relying on the convergence property (3.29) we have immediately 

n 

If = E( a K' h (*#-)V'(l#-)) - A J ' h (yJ-)) (u"' h ($-) - u'- h (yj-)) |u^(^+) - u^-)\ 

3 = 1 

so that 

n 

1^1 < c E E l« J,fc (^±) - «'fo±)l + l«" ,fc (i#±) - «"(»j±)I- 

3=1 ± 

Thus, in view of the local convergence at jump points (3.27a), for h small enough we obtain 

<Ce. (3.35) 

Relying on the simplifying assumption yj < y^ < yj +1 < yj +1 for all j, we can decompose T 2 h as 
follows: 

n „ 

T^ = Y^ / {a{u^ h ,u n ' h ) - / V*)) (u n ' h - u T ' h ) dV 1 * 

j=o •'(y l j>yj+i) 

- {aiu^u 11 ) - /'(u 1 )) (u 11 - u 1 ) dV 1 
~ E / M« V) - / V)) («" - «') ^ (3.36) 

n ,. 

+ E / (a^",^") - / V")) (</ J ' h - u*- h ) dV T ' h 

-•=o ■'[w + i.J/j+i) 



3 



-■ 1 2,1 ' 2,2 ' 1 2,3- 



We first consider T 2 h 2 : 



n „ 

r£ 2 = " E / {a{u I (y j +),u II (y j +)) - f'(n I (y j +))) {u II (y j +) - u I (y j +)) dV\y) 
j=o J (yj>Vj} 

n ,. 

+E / {(a(« 7 (»).«"(»))-/V(y))) («"(») -«'(»)) 

- (a(u J (y i +),u J7 (y i +)) - />'(%+))) (u"( yi +) - «'(%■+)) } dV'{y). 
Therefore, with (3.17), we obtain 

\T 2 h 2 \<cJ2\v I (yj+)-v I (yj+)\ 

3 

+ CV I (+^)( sup \u I (y)-u I (y j +)\+ sup \u H (y) - u II (y j +)\) . 
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Since yj — » yj, we have \ V I {yj+) — V I (yj+)\ — »■ 0, therefore for /i sufficiently small 

|T£ 2 |<Ce. (3.37) 



A similar argument for T^ 3 shows that 



Next consider the decomposition 



\T 2 h 3 \<Ce. (3.38) 



(a(u J '\ u H ' h ) - f(u^ h )) (u n ' h - u 1 *) dV z ' h - {a(u\u H ) - /'(u 1 )) (u 11 - u 1 ) dV 1 
= (a^'^u 11 ^) - /V'")) ~ « 7,h ) dV^ h - (aiu 1 ^ 11 ) - /V)) («" " «*) ^ 

+ (a(u J , u 11 ) - f'iu 1 )) (u 11 - u 1 ) dV T ' h - (a(ui , u 1 / ) - f'{u{)) {u 1 / - u{) dV T ' h 
+ (a(ui, u 1 /) - f\u{)) {u» - u{) dV 1 * - {a{u{,u{') - f(u{)) - u{) dV 1 
+ {a^lui 1 ) - f'{u{)) {ui 1 - u{) dV 1 - (afa'V) - /'(«')) («" " «') dV 1 , 

which, with obvious notation, yields a decomposition for T 2 ^ 1 

T 2 h 1 = Mi + M 2 h + M% + M 4 \ (3.39) 

Using (3.17) and the local convergence property (3.31), we obtain 



\M^\<C^2 \dV I,h \ ( sup \u I ' h -u I \+ sup \u J1 " - -u 



j=o^,Vi+i)' ' V (v?>Vi+0 (v?,w,-+i) 7 ( 3 - 40 ) 

< Ce. 



Similarly using (3.17) and (3.33) we obtain 

n n 

|M 2 h |<cV/ \dV z ' h \ ( sup |t/-uf| + sup |u 7/ -nf|) 

3=ohv^Vj+i) MvjNw+O (»*>Wi+i) 7 ( 3 - 41 ) 

< Ce. 



Dealing with M4 1 is similar: 



\M%\<C^2 SU P |^-^|+ sup |u /7 -u; 



1 -u{\+ sup |n /7 -nf 7 |) 
^Av?.wj+i)' ' V (i#,iy+i)' ' 7 ( 3 - 42 ) 



< Ce. 

Finally to treat M3 we observe that, since u{ and u 1 / are continuous functions on each interval 
(yj,yj + i) and since dV I,h is sequence of bounded measures converging weakly-star toward dV 1 , we 
have for all h sufficiently small 

\M%\ < e. (3.43) 

Combining (3.39)-(3.43) we get 

\ T ii\<Ce. (3.44) 
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Combining (3.36)-(3.38) and (3.44) we obtain 

|T 2 h | < Ce 

and thus with (3.34)-(3.35) 

|^2( r )| ^ C £ f° r an h sufficiently small. 
Since e is arbitrary, this completes the proof of (3.30). 

Step 3 : Consider now the term 

ni(r)= E (2m+TV{b h )^\a h (x-,T)-X h (x,T)\\u II ' h {x-,T)-u I ' h {x-,T)\. (3.45) 

On one hand, observe that 

TV{b h {T)) = TV{u^ h {T)) + TV{u n ' h (T)) — > V 1 (+oo, r) + F 7/ (+oo, r). (3.46) 

For all but countably many r the following holds. Extracting a subsequence if necessary we can 
always assume that for each j either (yj 1 , r) <G £(a h ) for all h, or else (y^ ,r) ^ ^(a' 1 ) for all /i. Then 
consider the following three sets: denote by J\ the set of indices j such that (y^r) G £(a ft ) and 
(yj,r) G £(ct). Let J2 the set of indices j such that (y^ ,r) £(a ft ) and (yj,r) G £(a). Finally J 3 is 
the set of indices j such that (y^,r) G £(0^) and (yj,r) ^ £(a). 
The local convergence property (3.27) implies 

■ ?eJl (3.47) 

— ► 5^l a K(fj—)» uJJ (fi-)) (%■+)) I -u 1 (vj-)\. 

jeJi 

(Indeed, given e > 0, choose finitely many jump points as in (3.31) and use (3.27) with e replaced 
with e\u I (z+) -u I (z+)\). 

On the other hand for indices in J2 or J3 we have 

E Hu I ' h {y^-),u II ^ h {y^-))-a{u I ' h {y l ;-),u I ' h {y^+))\ \u n - h {$-) -u J ' h (i#-)l — > (3.48) 
ieJ 2 uJ 3 

but 

E |aK(%-)V J (f;-)) -aK(y,-),n 7 (%•+)) | |^(y,-) - u'fo-)! = 0. (3.49) 

j€J 2 UJ 3 

Indeed, for each j G J2, yj is a Lax shock but y 1 - is not. Extracting a subsequence if necessary, it 
must be that the Lax inequalities are violated on the left or on the right side of y^ for all h. So it must 
be that, assuming that it is the case on the left side, o(u J {yj— ), u 11 (yj— )) — a(u 7 (yj —) , u 1 (yj+) > 
while o(n 7 (yj—),u H (y?— )) — o 1 (u I (y l j — ), u I (y^+) < for all /j. But the latter converges toward the 
former by the local uniform convergence, which proves that a(u 7 (yj — ) , u 1 1 (yj — )) — a^u 1 (yj — ), u 1 (yj+) =| 
0. 
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Combining (3.45)-(3.49) yields 

fif-*/ V q(r) \a(x-,r) - \(x,t)\ \u u (x,t) - u 1 (x,t)\ dr, (3.50) 



(x,T)eC(a) 

where g:=2m + V 7 (+oo) + V H (+oo). 

Step 4 : Continuity of the weighted norm. 

Fix some time t. Recall that the weight w(t) is defined based on the total variation functions 
V 11 and V 1 and on the function u IJ (t) — The weight w h {t) is defined based on the total 

variation functions V H ' h and V I,h and on the function u n ' h (t) — u I,h (t). On the other hand, 
U n,h _ u i,h( t ) u 11 - u 1 ^), V u ' h -f F 7/ and V J ' fc -» y J . Therefore we have 

w(x, t) = w(x, t) whenever u II (x, t) — i/(x, t) ^ 0. (3.51) 

Combining (3.51) and the L 1 convergence u H ' h — u I,h (t) — > n /J — u J (t), we have 

\\u H (t) - ^(t)|U (t) = lim \\u n > h (t) - u^ h {t)\\ wh{t) . (3.52) 



Step 5 : The left-hand side of (3.16) is bounded above by (3.21). 
First of all, the inequality 

/* / (a(« J , u 11 ) - fin 1 )) (u 11 - u 1 ) dV\y, r)dr 
Jo Jr 

<ff (a(« / ,« /J )-/ , (« J ))(« /J -« / )^ / (j/,r)dr 
Jo Jr 



(3.53) 



is a direct consequence of (3.19) and the definition of the nonconservative product in (3.15). 

On the other hand, by the definition of the weighted norm and because of (3.20), similarly to 
(3.51) we have the inequality 

w(x,t) < w(x,t) whenever u II (x,t) — u z (x,t) ^ 0. (3.54) 

Hence, by (3.53) and (3.54) the left-hand side of (3.16) is bounded above by (3.21). This completes 
the proof of Theorem 3.5. □ 

4. Generalized Characteristics and Maximum Principle 

We now return to the setting in Section 2 and aim at extending the analysis therein to arbitrary 
functions of bounded variation. For exact solutions of the hyperbolic equation 

d t ^ + d x ( ai p) =0, (4.i) 

we will establish a maximum principle: Any solution of (4.1) remains non-negative for all times if 
it is so initially. For a more precise (local) statement, our proof will make use of Dafermos-Filippov 
theory of generalized characteristics. 
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Our main assumption throughout this section is the following: 

E 

There exists a constant E such that d x a < —. (4.2) 

This is nothing but a generalization of the well-known Oleinik's entropy inequality. To motivate 
(4.2), let us recall the following result. 

Let / be a strictly convex function and u be an entropy solution (with bounded variation for all 
times) of the conservation law 

d t u + d x f{u) =0, u(x,t)eR. (4.3) 
Then is is known that there exists a constant C = C(u) such that 

d x u<j. (4.4) 

Lemma 4.1. If u 1 and u 11 are two entropy solutions of the conservation law (4.3), then the aver- 
aging speed 

a = a(u I ,u») := f f ] ■ (^) 

u 11 — u 1 

satisfies our assumption (4.2), with E = sup/" (C(u L ) + C(u II ))/2. 

Proof. Let us fix some time t > 0. On each Borel set consisting of points of continuity of both u 1 
and u 11 , the following holds: 

d x a = d x [ f'(9u I + (l~9)u n )d9 
Jo 

= f f" (eu 1 + (i - e) u 11 ) (e d x u* + (i-o) d x u H ) 

Jo 1 1 

Ciu^ + Ciu 11 ) 



de 



< sup/ 



2t 



On the other hand, at a point x where one of u 1 or u 11 is discontinuous, we have with an obvious 
notation 



a_|_ — a_ 



/ f(9u I + + (l-9)u I + I )de- [ f'(9u I _ + (1 - 9)u^) d9 < 0, 
Jo Jo 



since /' is an increasing function and (for instance by (4.4)) both u 1 and u 11 satisfy < u{_ and 
u 1 ^ 1 <uU- □ 

By definition, a generalized characteristic y = y(t) associated with the coefficient a must satisfy 
for almost every t (in its domain of definition) 

a+(y(t),t)<y'(t) <a-(y(t),t). (4.6) 

According to Filippov's theory of differential equations [12], through each point (x, t) there pass a 
maximal and a minimal generalized characteristic. 
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Definition 4.2. A generalized characteristic is said to be genuine iff for almost every t it satisfies 

y'{t) G {a-(y(t),t),a + (y(t),t)}. (4.7) 

Proposition 4.3. Any minimal backward generalized characteristic is genuine and for almost every 
t satisfies 

y'(t) = a.(y(t),t). (4.8) 

Similarly, for a maximal backward generalized characteristic we have y'(t) = a + (y(t),t). 

Proof. Here we only rely on the following consequence of (4.2): a + < a_ at each discontinuity point 
of the function a. Geometrically, this condition prevents the existence of rarefaction-shocks in a. 
On the other hand, rarefaction centers (also prevented by (4.2) for t > 0) could still be allowed for 
the present purpose. 

Consider (x,t) G (— oo, +oo) x (0, oo), and let y(t) := y(t;x,t) be the minimal backward charac- 
teristic through (x,t). We prove that it is genuine on its domain (s,t\. We proceed as in [8] and 
assume by contradiction that there is a measurable set J, J C (s, t\ of positive Lebesgue measure, 
and e > such that 

a-(y(t),t)-y'(t)>2e, t G J. (4.9) 

For each t G J there exists 5(t) > with the property 

a+(x,t)>a-(y(t),t)-e, x G (y(t) - S(t), y(t)). (4.10) 

Finally, there is a subset I C J with /i* (I) > (here fi* denotes the outer measure) and 5 > such 
that S(t) > 5 for i G /. 

Let t be a density point of /, with respect to fi* . Thus there exists f, < f < i — r, so that 

fi*(ln[T,r + r]) 2\a\+e n ^_ 

— — > -r-T , < r < r, (4.11) 

r 2|a| + 2e' ~ y J 

where 

a := inf{o + (x, t) — a-(y(t),t) : s <t <t, y(t) — 5<x< y(t)}. 

Now take a point y G (y(r) — 5, y(r)) with the property y > y(r) — \er, and consider a forward 
characteristic z(-) through (y,r). We first observe that 

z(t) < y(t), t > r, 

since y(t) is the minimal backward characteristic through (x,t). 
In addition, we have 

z(t) > y(t)-S, t >G [r, T + f}. 

Indeed, suppose by contradiction that for some r G (0, f], z(t) > y(t) — 5 for t >G [r, r + r), but 
z(t + r) = y(r + r) — 5. Then 

z'(t)dt - j/(t) - y + 6 

> J T+ \z'{t)-y'{t))dt 

= I (z'(t) - a_(y(t),t) + a.(y(t),t) - y'(t))dt 

Jln[T,T+r] 

+ [ (z'(t) - a-(y(t),t) + a-(y(t),t) - y'(t))dt 

J\T,T+r]\I 



r]\ 

> efi*(lD [r,r + r]) +a(r- n*(lD [r,r + r])) > 0, 
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by (4.9)-(4.11), which leads to a contradiction. In the same way one obtains 



z'(t)dt - y(t) - J y'(t)dt 
> £H*(ln[T,T + f]) + a(f - n [r,r + r])) - \ef > 0, 



which gives another contradiction. For the maximal backward characteristic the proof is similar. 
□ 

Proposition 4.4. Forward characteristics leaving from some (x,t) are unique when t > 0. 

Proof. Suppose there were two forward characteristics y(-) and z(-) through (x,t) with y(r) < z(t) 
for some t > t. By (4.2) we have 

z'(t) - y'(r) < a-(z(t),t) - a + (y(t),t) < C- t {z(r) - y(r)). (4.12) 

Integrating (4.12) from i to r one gets z(t) — y(r) = 0, which gives a contradiction. □ 
Theorem 4.5. Let tp = ip(x,t) be a solution of (4.1) such that on some interval [£o,Co] we have 

#r,0)>0, x€[&,Co]- (4.13) 

Let £ = ^(t) be any forward generalized characteristic leaving from (£o>0) ; and ( = ((t) be any 
forward generalized characteristic leaving from (£o, 0). 
Then we have for all t > 

*P(x,t)>0, xe{m,C(t)). (4.14) 

Note that it may happen that £(i) = ((t) for t large enough. 

Proof. Observe that the two characteristics cannot cross and fix any time t > such that £(t) < CW- 
Fix also any two points such that £(t) < y < z < C(t)- Let y(t) and z{t) be the maximal and minimal 
backward characteristics emanating from y and z, respectively. These characteristics can not leave 
the region limited by £(t) and ((t). 

Integrating (4.1) in the domain bounded by the characteristics y(t) and z(t), and using that these 
characteristics are genuine, so that the flux terms along the vertical boundaries vanish identically, 
we arrive at 

,z ^(0) 

/ ^(x,t)dx= / tp(x,0)dx>0. (4.15) 

Jy Jy(0) 

The last inequality is due to the fact that if)(., 0) > and the inequalities = < 2/(0) < z (0) < 
C(0) = Co- Since y and z are arbitrary, we obtain (4.14). □ 

5. A Sharp L 1 Estimate for Hyperbolic Linear Equations 

Based on the maximum principle established in Section 4, we now derive a sharp estimate for the 
weighted norm introduced in Section 2. We restrict attention again to the situation where u 1 and 
u 11 are two entropy solutions of the conservation law (4.3) and a is the averaging speed given in 
(4.5). We define a weight by analogy with what was done in Section 2 in the special case of piecewise 
constant solutions. 
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Given a solution -0 of the equation (4.1), we introduce weighted L 1 norm in the following way. 
Set 

V I (x,t)=TVf 00 (u I (t)), V II (x,t)=TV* 00 (u II (t)) (5.1) 

and fix some parameter m > 0. Then consider the weight-function defined, for each t > and each 
point of continuity x for u 1 \t) and u 11 (t), by 

{m + V T (oo, t) - V T (x, t) + V n (x, t) if i/>(x, t) > 0, 
(5.2) 
m + V 7 (x, t) + V 11 (oo, t) - y /7 (x, t) if i/j(x, t) < 0. 

It is immediate to see that 

m<w(x,t)<m + TV(u I (t))+TV(u II (t)), x G R. (5.3) 
Finally the weighted norm on the solutions ip of (4.1) is defined by 

\\^(t)\\w(t) ■■= / \ip(x,t)\w(x,t)dx. 

JR 

Note that the weight depends on the fixed solutions u 1 and u 11 , but also on the solution V 1 - 
Our sharp estimate will involve the nonconservative product 

^(t) = (a-f'(u I (t)))mdV I (t) 

defined for all almost every t > by 

(1) If B is a Borel set included in the set of continuity points of u 1 (t) then 

^{t){B) = [ (a(t) - / V(t))) dV'(t), (5.4a) 

where the integral is defined in a classical sense; 

(2) If x is a point of jump of ii J (t), then 

/4(*)({*}) = (a(x-,t) - A J (x,t)) V(x-,t) |n J (x+,t) - u J (x-,t)|. (5.4b) 

Here A / (x, t) is a the shock speed of the discontinuity in u 1 located at (x,t). The measure ^ (t) 
is defined similarly. Regarding the expression (5.4b), it is worth noting that if (x,t) is a point of 
approximate jump of u 1 and then the jump relation for the equation (4.1) reads 

(a(x-,t) - A 7 (x, t)) ij)(x-,t) = (a(x+, t) - A'(x, t)) ^(x+, t). (5.5) 

In the same way we define 

4 I (t) = (f'(u II (t))-a)mdV II (t). 

We now prove: 
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Theorem 5.1. Let u 1 and u 11 be two entropy solutions of (1.1) such that u 11 — u 1 admits finitely 
many changes of sign. Let ip be any solution of bounded variation of the hyperbolic equation (4.1) 
satisfying the constrain 

ip (u 11 - u 1 ) > 0. (5.6) 



Then for all < s < t 

W {t)+ [ (^m + TV(a))\a.(x,T)-X(x,T)\\^_(x,T)\dr 

Js (x,r)eC(a) 

+ f I (a(T)-f\u I (T)))iP(T)dV I (r)dT+ f [ (a(T)-f'(u H (T)))^T)dV n (T)dT 

Js Jr Js JR 

< ll^(a)IU( fl )- 

(5.7) 

The assumption (5.6) is clearly satisfied with the choice ip = u 11 — u 1 . Therefore our previous 
result in Theorem 3.5 (derived via a completely different proof) can be regarded as a corollary of 
Theorem 5.1. 

It is interesting to observe that, when u 11 = u 1 , the weight (5.2) becomes constant, and therefore 
(5.7) reduces to the L 1 estimate. 

V(R)+ [ (2m + TV(a)) |a_(x,r) - X(x,t)\ ^.(x,r)\dr 

< \Ms)\\ LHR) . 

Also, note that under the assumption (5.6) /U^(i) and (t) are positive except at points (x,t) € 
C(a) U lZ(a). However, these negative terms are offset in (5.7) by the positve terms under the first 
integral. 

Proof. Fix any positive time t. By assumption we have finitely many points — oo = y < y\ < 
. . . < y n < y n +i = +°° such that, on each interval (yj,?/j+i), we have ip(t) > when i is odd and 
ip(t) > when i is even. For every i = 1, • ■ ■ ,n, consider the (unique by Proposition 4.4) forward 
characteristic yi(-) associated with the coefficient a and issuing from the initial point (yi,t). 

We will focus attention on some interval (j/i, y^+i ) with i odd, say, and with — oo < < yi+i < 
+oo. Except when specified differently, all of the characteristics to be considered from now on 
are associated with the solution u 11 . For definiteness we will first study the case that the forward 
characteristic Xo(-) (associated with u 11 and) issuing from the point (yi,t) is located on the right-side 
of the curve yi, that is, 

2/;(r) < xo(r), t<r<t + 5 

for some 5 > sufficiently small. 

Fix some (sufficiently small) e > and denote by yi < Z\ < . . . < < yi+i the points where u 1 
has a jump larger or equal to e, that is, 

u I _ I (z I ,t)-u I + I (z I ,t)>e, I = 1,...,N. (5.8) 

For each / = 1, . . . , N, consider also the forward characteristic Xi{') issuing from the point (zj, t). 
For definiteness, we will also assume that the forward characteristic xn+i{-) issuing from (y i+1 ,t) 
satisfies 

Xn+i(t) < y i+ i(r), t<r<t + 5 
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for some 5 > sufficiently small. 

Next, let us select a time s > t with s — t so small that the following properties hold: 

(a) No intersection among the characteristics Vi,Xo,Xi-> ■ ■ ■ ,XN,XN+i,Vi+i occurs in the time 
interval [t, s]. 

(b) For 1=1,... ,N,let ("/(•) and £/(•) be the minimal and the maximal backward characteristics 
emanating from the point (x/(s),s). Then the total variation of u H (-,t) over the intervals 
(Q(t),zi) and (zi,£i(t)) should not exceed jj. 

(c) Let Co(') be the minimal backward characteristic emanating from (yi(s),s) and £o( - ) be the 
maximal backward characteristic emanating from (xo{s),s). Then the total variation of 
u n (-,t) over the intervals (yi,^o(t)) and (Co(t),Vi) should not exceed e. 

(d) Let Cat+i(-) be the minimal backward characteristic emanating from the point (xn+i(s), s) 
and £jv+i(-) be the maximal backward characteristic emanating from (y i+1 (s), s). Then the 
total variation of u H (-,t) over the intervals (£jv + i(i), Vi+i) an d (yi+ii^N+i(t)) should not 
exceed e. 

For 1 = 0,... ,N, and some integer k to be fixed later, consider a mesh of the form 

Xi(s) = x? < x) < . . . < x k j < x k z +1 = xi+i(s). (5.9) 

For / = 0, . . . ,N and j = 1, . . . , k, consider also the maximal backward characteristic ema ~ 
nating from the point (x J j,s) and identify its intercept zj = (,j(t) by the horizontal line at time t. 
Finally set also 

z° =y l , z k N +1 = y i+1 , z k +l = z° T = z u I=1,...,N. 

To start the proof, we integrate the equation (4.1) satisfied by the function ip, successively in 
each domain limited by the characteristics introduced above. Applying Green's theorem, we arrive 
at the following five formulas: 

(i) Integrating (4.1) on the region 

{{x,r)/t < t < s, yi{r) < x < xq{t)} 
and multiplying by V II (yi,t) one gets 

-Xo(s) 



i{j(x, S )V n ( yi ,t)dx+ / (y' t ~a + )i; + (y l (T),T)V II (y l ,t)dT 

+ [( a -- A o) ^- (Xo(r), r) V H ( yi , t) dr = 0. 
(ii) Integrating (4.1) on each of the regions 

{(x,r)/t<r<s, e^KxK^ir)} 
for 1 = 0,... , N and j = 1, . . . ,k, and then multiplying by V 11 (zj+, t), one gets 



(5.10i) 



3+1 3+1 

c i . r z i 

ip(x, s) V n (z 3 j+, t)dx- / ip(x, t) V II {z 3 I +, t) dx 



s 

J ay 

+ ^V,-g+)V'+(#(t),t)V"(4+^ 

(5.10ii) 
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(in) Integrating (4.1) on each of the regions 

{(x,r)/t < t < s, xi (r) <x<£(t)} 
for / = 0, . . . ,N, and multiplying by V II (zj+, t) one gets 



r x i r z i 

/ ip(x,s)V II (z I +,t)dx- / '4}{x,t)V II {z I +,t)dx 

Jxi(s) Jzi 

+ J (\i~a + )i; + (xi(T),T)V II (z I +,t)dT + J (a. - \ 1 I )iP-(ei(T),T)V n \z I+ ,t)dr = 0. 

(5.10iii) 

(iv) Integrating (4.1) on the regions 

{(x,r)/t < t < s, $(r) < x < Xi+i(t)} 
for I = 0, . . . ,N, and multiplying by V H (z k +, t) one gets 
rxi+i{s) rzi+i 

/ i){x,s)V II {z k I +,t)dx - / ^{x,t)V II {z k I +,t)dx 

+ | S (A^-a + )V+(ef(r),r) V H (z k +,t) dr + J\a_ - A /+1 ) V-(x/+i(r),r) V 7/ (zf +, i) dr = 0. 

(5.10iv) 

(v) Finally integrating (4.1) on the last region 

{(x,t) /t < t < s, Xn+i(t) < x < y i+ i(r)} 
and multiplying by V n (yi + i, t) one gets 

/S/i+i(s) /.s 
ip(x,s) V H (y i+1 ,t) dx + / (Aat+i - a+) ^ + (xjv+i(r),r) V" JJ (j/ i+ i, t) <Jt 
~*+iW -/t (5.10v) 

+ j {a--y[ +l )^-{y t+l {T),T)V n {y t+l ,t)dT = V. 

Next, summing all of the formulas (5.10) leads us to the general identity: 

rxo(s) N k f x 3 j +1 



'vi( s ) l=0j=0" x / 

rlK+i(s) W fc + 1 



+ f + ^(x,s)V I (y i+1 ,t)dx-J2J2 I ' ^(x,t)V n (zi+,t)dx 

N k s 

= EE / [V II (zi+,t)-V II (zr 1 +,t)}(X^-a.)^(e i (r),r)dT 

1=0 j=l Jt 

"E / [V H (z I+ ,t) - V n (z k _i+,t)}^ - a.)^( X i(r),r)dr 

- J [V u ( yi +, t) - V n ( yi ,t)} (A - a_) V- (xo(r), r) dr 

- [F /7 (y l+1 +,t) - V 7/ (4,t)] (A w+1 - a_) V-(Xiv+i(r),r) 

- jf (y- - a+) V+(2/i(r), r) V 7 %, t) dr - J (a_ - ^_(y l+1 (r), r) V /7 (y l+ i, i) dr. 
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To estimate the right-hand side of (5.11), we recall that the solution u 1 of a scalar conservation laws 
satisfies 

V n (y t ,t) > V H ( X o(s),s), V n (zj+,t) > V n {x{+,s), 

for I = 0, . . . , N and j = 0, . . . , k. Hence, choosing the difference x\ +1 — x\ in (5.9) sufficiently 
small and since the function V H (-,t) is nondecreasing, we conclude that the left-hand side of (5.11) 
can be bounded from below, as follows: 

ryi+i{s) ryi+i(t) 
L.H.S. > / ^(x,s)V n (x,s)dx-(s-t)e- ip(x, t)V n (x, t)dx. (5.12) 

■'yi(s) Jyi(t) 

Estimating the right-hand side of (5.11) is more involved. First note that each term arising in 
the left-hand side of (5.11) is non-positive. This follows from our condition (5.6). Indeed, consider a 
point (x, s) of approximate jump or approximate continuity of u 1 , u 11 and ip. If all of these functions 
are continuous, the result is trivial. Call A the discontinuity speed. Based on the jump relation (5.5), 
we see that either ^_ (A — a_) = ip + (A — a+) = 0, or else all of the terms ip-, A — a_, ip + , and 
A — o_|_ are distinct from zero. 

Suppose first that (x, s) is a point in the interior of the region limited by the two curves yj(.) and 
yi_l_i(.). In the latter case, since tp > in the region under consideration, we deduce that ip- > 
and ip + > 0, while the terms A — a_ and A — a + are either both negative or both positive. Actually, 
in view of the sign condition (5.6), we have u± — u± > and, therefore, A — a± > as follows from 
(3.10) (here we are dealing with a jump of u 11 ). 

Consider next a point of the boundary t/i, for instance. So we now have ip- < and ip+ > 0, while 
the terms A — a_ and A — a + opposite sign. Since no rarefaction-shock can arise, the discontinuity 
must be a Lax shock and so A — a_ < and A — a + > 0. Again the corresponding term in (5.11) 
has a favorable sign. (Observe that the condition (5.6) was not used in this second case.) 

Then, for all / = 0, . . . ,N and j = 1, . . . ,k, let #}(•) be the (maximal, for definiteness) backward 
characteristic associated with u 1 and issuing from the point (£j(t),t). Denote also by 9(z j i ;t) its 
intercept with the horizontal line at time t. Setting 

~, , , _ /(^(x,t))-/(^(fl(x,T),t)) 

l ' ' J ' u II (x,t)-u I (6{x;T),t) 

and using that the solution u 1 remains constant along the characteristic Oj(-), we obtain 

(\>-a-)(e i {T)) = \i{4)-a(4 1 t-,T). (5.13) 

Then consider the (maximum, for definiteness) backward characteristic yj(-) associated with a and 
issuing from the point (£j(t),t). By integrating tp along the characteristic y\{ ) and using the 
inequality (4.4), we arrive at a lower bound for ip 

m(T),T)>^{y]{t),t)(^j , t<r<s. (5.14) 

Upon choosing x{ +1 — x\ in (5.9) so small that the oscillation of V^{-) over each interval (zj — zj +1 ) 
does not exceed e and recalling the standard estimates on Stieltjes integrals we deduce from (5.11)- 
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(5.13) that 

N k s 

EE / [V II (4+,t)-V II (4- 1 +,t)} ((Aj-a_)V-(ej(r),T))dr 

N k fS / f\ E 

^EE/ [K"(4,*)-K"(^r\0](Ai(4)-a(4,t,r))^(i),t) {-) dr 

1=0 3 = 1 Jt W 

(Aj(x) - a(x,t,T))i/)(x,t) dV c n {x,t) - c{y i+1 - yi )ej I -J dr. 



(5.15) 



t \Jy 



We now combine (5.10), (5.11) and (5.15), divide the resulting inequality by s — t, and let s\t, 
e — > 0, obtaining the following inequality: 



rVi+i(t) rVi+i 

/ ^(x,t)V H (x,t)dx < - / (A J 7 (x) - a(x, t)) ^(x, t) dV^(x, t) 

Jyi(t) Jyi 



~dt 

- E («_(ar,t)-<(a:,t))(A J -a_)(x,t)V-(ar,t) 
(i,t)ej(ti") 

- -u^( yi ,t)) (X 1 - a-)(yi,t)ip-( yi ,t) 

- (uL 1 (y i+1 , t) - u+ (y i+1 , t)) (X 1 - a_) (y i+ i, t) V- (Vi+i, t) 
-{y' i -a + )^+{y h t) V n { yi ,t) 

- (a- -y' i+1 )ip-(y i+1 ,t) V H (y l+1 ,t) 

(5.16) 

The third and fourth terms in the right-hand side of (5.16) are due to the fact that xo and Xn+i he 
inside the region limited by yi and j/j+i- 

We can next focus on the intervals (y^, t/j+i) with i even. Based on a completely symmetric 
argument and using now the weight m + V 11 (oo, t) — V 7/ (-, i) instead of V IJ (•, t), we obtain 



< 



j+ rVi+iii) . N 

' / (-^(x,t)) (m + y 7/ (oo,t) -y /7 (x,t))cix 

A/«(t) V 7 

(A^(x) - a(x, t)) (-^(x, t))^/ / (x, t) 

+ E (*>') ~ u + (*>*)) (A 7 -a_)(-^_)(x,t) 

OMjeW) (5.17) 
+ {u-(yi,t)-u^(yi,t)) (A / -a_)(-^_)(y,,t) 

+ («i J (y i+ i,t) -< J (y i+ i,t)) {X 1 - a-)(-tp-)(y i+ i,t) 
-(y[-a + )(-^ + )(y i ,t){m + V II (oo,t)-V II (y i ,t)) 
-(a_-y< +1 )(-V-)(l/i + i,t)(m + (oo,i) - V% +1 ,i)) . 
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By summation over i = 1, . . . , re in (5.16) for i odd and in (5.17) for i even respectively, we obtain 



dt 



[V(ar, t)] V /7 (x, t) + [-ip(x, t)] + (m + V n {oo, t) - V n {x, t)) dx 



<- ^2 ( m + V n (oo,t)^ \\{x,t) - a-(x,t)\ \ip-(x,t)\ 

{u-(x,t) -u^(x,t)) (A / (x,t) -a-(x,t)) i[>-(x,t) 

(x,t)ej(u") 

(f'(u n (y,t)) - a(y,t)) ^(y,t) dV. 11 (y,t), 



(5.18) 



R 

where the superscript + denotes the positive part of the functions -0 and — -0 respectively. 
Consider now the case where 

Xo(t) < ^(r), t<r<t + 5, 

and 

Vi+i{r) < Xn+i(t), t<T<t + d. 
Assume that there exists a time f > t such that 

Xo(r) < yi(f), Vi+i{f) < XN+i(f) 

(otherwise the curves of the two pairs will coincide, and we can reduce to the previous case). Let 
now £o( - ) be the maximal backward characteristic emanating from (yi(f),f), and (jv+i(-) be the 
minimal backward characteristic emanating from the point (yi + i(f ), f ). Since characteristics cannot 
cross, we have that 

Vi{t) < £o(i), Civ+i(i) < y i+ i(t), 
Then, by finite propagation speed, there exists a time s > t such that 

Vi{r) < £o(t), (n+i(t) < y i+ i(r), t<r < s, 

Vi{s) = £o(s), (n+i(s) = y i+ i(s). 

Instead of properties (c), (d), we will require that s satisfies the following: 

(c') Let Co(') be the minimal backward characteristic emanating from (xo(s),s). Then the total 

variation of u /7 (-,i) over the intervals (yi,£o(^)) and (Co(£)>2/i) should not exceed e. 
(d') Let £/v+i(-) be the maximal backward characteristic emanating from (xn+i(s), s). Then the 
total variation of u H (-,t) over the intervals ((N+i{t),yi+i) and (y i+1 ,S t N + i(t)) should not 
exceed e. 

From then on we can proceed as before. Finally we write the inequality in (5.18) exchanging the 
roles of u 1 and u 11 , and combining it with (5.17) we arrive exactly at the desired inequality (5.7) 
and the proof of Theorem 5.1 is completed. □ 
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